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O I We study the one-loop effective action for 4D conformally invariant 

Q>^ I scalars on the spherically symmetric background. The main part of effective 

action is given by integration of 4D conformal anomaly. This effective action 
(in large approximation and partial curvature expansion) is applied to in- 
vestigate the quantum evolution of Schwarzschild-de Sitter (SdS) black holes 
r-| ! of maximal mass. We find that the effect (recently discovered by Bousso 

and Hawking for minimal scalars and another approximate effective ac- 
is-H " tion) of anti-evaporation of nearly maximal SdS (Nariai) black holes takes 

■ also place in the model under consideration. Careful treatment of quan- 



tum corrections and perturbations modes of Nariai black hole is given being 
quite complicated. It is shown that exists also perturbation where black 
hole radius shrinks, i.e. black hole evaporates. We point out that our result 
holds for wide class of models including conformal scalars, spinors and vec- 
tors. Hence, anti-evaporation of SdS black holes is rather general effect which 
should be taken into account in quantum gravity considerations. 
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1 Introduction. 



The gravitational black holes radiation [|T| may be considered as the one of 
the most brilliant manifestations of quantum gravity. Moreover, one can 
expect that other quantum gravity strong effects should be also searched in 
black holes. Recently, one such new effect (anti-evaporation of black holes) 
has been found in ref. 0. The authors of above work used large and s- 
waves approximations in their calculation of effective action for 4D minimal 
scalar. This effective action has been used to study the quantum stability of 
the Schwarzschild-de Sitter black holes of maximal mass (Nariai solution). 
It has been shown that there is the perturbation mode where the black hole 
size increases and the black hole grows back to the maximal radius (anti- 
evaporation). 

The interesting question is if that effect is quite common or it results from 
the specific model and (or) approximation used. To answer this question we 
study another model (of 4D conformal scalars) and find the effective action 
in another approximation (large and expansion on the curvature). The 
analysis of quantum corrected equations of motion is much more complicated 
than in Bousso-Hawking model. Nevertheless, the final result is quilitatively 
the same- we confirm the possibility of black holes anti-evaporation. 

The contents of the paper are following. In the next section 4D anomaly 
induced effective action for conformal scalar is evaluated. We work in large 

approximation and calculate the piece of the effective action which is not 
anomaly induced as expansion on the curvature near spherically symmet- 
ric background. As a result we obtain effectively 2D gravitational system. 
Equations of motion are derived in the section 3. Section 4 is devoted to the 
analysis of quantum perturbations for Nariai black hole. We demonstrate 
the possibility of anti-evaporation of such objects. In the last section we give 
some comments and explain that our result is of general relevancy for a class 
of quantum conformal matter models. 

2 Effective action for conformal scalars 

In the present section we derive the effective action for conformally invariant 
scalars (for a general review of effective action in curved space, see[|]). Let 
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us start from Einstein gravity with conformal scalars Xi'- 
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The convenient choice for the spherically symmetric spacetime is the follow- 
ing: 

ds'^ = g^^dx^'dx" + /(0)c/fi (2) 

where ^,1^ = 0, 1, g^^, and depend only from x^. 

Reducing 4D action (|1|) in accordance with the metric (0) and working in 
s-wave sector (i.e., Xi = Xi{x^i^^))y we get 



S_ 

An 



N 



d X\f^^ 



f{<P){R-2A) + Rn + 2{V^p){VJ 
f{ct>)V"x^^o.X^ + \xl{Rf{^) + Rn) 





(3) 



where factor An appears as volume of S"^ space and i?n = 2 (or i?^ = ^ 
where = 1). If we take i?n = ^, then "1/53 = Anp^. 

Let us start the calculation of effective action due to scalars on the back- 
ground (^. 4D scalar in Eq.(|l]) is conformally invariant. Let us rewrite the 
metric (|^) as following: 

ds^ = f{(P) [f-\(f))gf,^dx^'dx'' + dn] . (4) 
In the calculation of effective action, we present effective action as following: 

r = r,„, + r[i,^7it^] (5) 

where Tind = T[f, g^^J] — r[l, gj^J] is conformal anomaly induced action which 
is quite well-known gj^J is metric (^) without multiplier in front of it, i.e., 
g'^J corresponds to 



ds' 



r\(f))g^^dx''dx'' + dn 



(6) 
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The conformal anomaly for 4D scalar s is well-known one: 



T = b{F + -nR] + b'G + b"nR 



(7) 



where b 



N 



b' 



N 



b" = but in principle, b" may be changed 



120(47r)2 ' 360{47r)2 ' 

by the finite renormalization of local counterterm, F is the square of Weyl 
tensor, G is Gauss-Bonnet invariant. 

Conformal anomaly induced effective action Find may be written as follows 



W 



a 



bjd^x^gFa 
+b' J d!^x^J—g^a 

+ (g-?m).} 

+ + j d'x^lR-ma (8) 

where a = |ln/(^), and a-independent terms are dropped. All 4- 
dimensional quantities (curvatures, covariant derivatives) in Eq.(|D should 
be calculated on the metric (|^). (We did not write subscript (4) for them.) 
Note that after calculation of (^ on the metric (^, we will get effectively 
two-dimensional gravitational theory. 

The next step is to calculate second term in the right hand side of Eq.(H), 
i.e., conformally invariant part of EA. We work in s-wave approximation 
where scalars depend only from first two coordinates. So we have to make 
reduction for scalar action on the metric (|^). We present now the metric (|^) 
as follows: 

ds"^ = g^^dx^'dx" + dVL (9) 

where g^y = f~\(j))g^,y. 

On the space with the background (^, having in mind s-wave approxi- 
mation and reduction, we get the following two dimensional classical action 
for scalars: 



S 
4tt 



2 E y d'x^g^"x^^^x^ + gX- {R + 



(10) 
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where all quantities in (|^) (i.e., R, □ and V^) should be calculated on two- 
dimensional metric g^^. As one can easily see second term in (|l^) breaks the 
conformal invariance. Denote '^"t^" = V. 

D 

It is clear that it is impossible to calculate the finite part of effective 
action corresponding to quantum system (^) in closed form. So, we will 
develop expansion on the "potential term" V. Zero term of this expansion 
(i.e., when we put = 0) may be found in closed form as it is standard 
Polyakov anomaly induced action 0: 

The first order term on the potential V may be found using zeta- 
regularization technique (see, ^ for a review). Actually, this term represents 
Coleman- Weinberg type term: 

r,., = ^/Ay=§(v-Vln^) (12) 

where /i^ is an arbitary mass parameter. Note that one can calculate the 
next terms of the expansion on V using different techniques. If we suppose 
that curvature R is quickly changing than next leading term of expansion is 
of the form ( pH] ) with another (significantlly smaller) coefficient than in (|Tl|). 
Finally effective action in our model will take the following form: 

r = Tind - r(o) - r(y) (13) 

where Find is given by Eq. and all quantities in Eq. (^ should be calculated 
on metric (H), r(o) is given by two-dimensional expression ([l^) and T(^v) is 
given by Eq . (|12D . The quantum effective action (|13D should be added to 
classical action (|^). That is combined action iS + F which will define the 
quantum dynamics of the system under discussion. 

We now consider to solve the equations of motion obtained from the above 
effective Lagrangians 5" -|- F. In the following, we use Qf^u and cr as a set of 
independent variables and we write g^u as g^u if there is not confusion. First 
we should note that the definitions of the Gauss-Bonnet invariant G and the 
Weyl tensor C^^a/s are given by 

G = R^,a(3R'""'^ - ^RapR""^ + R^ (14) 
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+ ^R{9f,a9ul3 - 9^lf39ua) ■ (15) 

Then the square of the Weyl tensor has the following form: 

Cfj.ualsC^""'^ = Rfiua/sR^""'^ " '^Ra/sR'^^ + -R^ (16) 

Then Find (W in Eq.(||)) is rewritten after the reduction to 2 dimensions as 
follows: 

+^(V'^/?(2))V,) a 
2i?^i?(^) - ^□i?^^)) a} 

X + i?n - 6na - 6V^aV,,a)' - (i?^^) + . (17) 

Here Rq_ = 2 is scalar curvature of with the unit radius. The suffix "(2)" 
expresses the quantitiy in 2 dimensions but we abbreviate it if there is no any 
confusion. We also note that in two dimensions the Riemann tensor R^uap 
and R^u are expressed via the scalar curvature R and the metric tensor g^y 
as follows 

R^lV(Tp = - {.9pa9vp — 9pp9ua) R , Rpv = -^9fiyR ■ (18) 
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3 Equations of motion 



Let us derive the equations of motion with account of quantum corrections 
from above effective action. In the following, we work in the conformal gauge 

^±T = -^e2^ ^±± = (19) 

after considering the variation of the effective action F + with respect to 
g^y and a. Note that the tensor g^^, under consideration is the product of 
the original metric tensor in and the a-function e~^'^, the equations given 
by the variations of g^^^ are the combinations of the equations given by the 
variation of the original metric and a-equation. 

Since g^^ and g vanish in the conformal gauge (p!9|) , we can drop 
second or higher order terms in g^^ and g and the scalar curvature R and 
de Alembertian □ have the following forms: 

8e"2^a+9_p 



R 

Q 



4g 

g- (-dlp-2{d^p)'' 



-^d_g - -d-g 



d-p 



^d+g+^d+p , 



+^/++ [-dlp-2{d^pr 
□ = -Ae-^Pd+d.+e-^Pd+{g++e^''d+)+e-^''d_{g—e^Pd_) 
In the conformal gauge (|19|), the action Vind ( |T^ ) is reduced into 



(20) 



J- ind 

47r 



32 



-4p 



16 



P 



a 



+b' J d^xe^p{a [l6e-^Pd+d-{e-^Pd+d- 



cr 



+e 



-2p 



^-^e-''>d^d^p+f]d^d.a 



2" [d+ie~^fd+d^p)d^a + d^{e~^Pd+d^p)d+a 

+ { 16e-^Pd+d^p + ye-2pa+a_(e-2^a+<9_p)^ a| 
2 



X [2{d+d^a + d+ad^a)e-^P{8e-^Pd+d_p + 2) 
+24e-'^f{d+d_a + d+ad_a) 



(21) 
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and the variations of Find with respect to g^^ are given by 



b' 



+-e^Pd±ad±{{8e-^Pd+d^p + 2)a} 
3 

+-e^Pad±i8e-^Pd+d^p + 2)d±a 
3 



+-e^^d±{8e-^Pd+d.p + 2)d±a 



- |6" + -{b + b')j [4e^Pd±ad±{8e-^''d+d-p) 

-A{d±afd+d-p + l2e^''d±ad±{e-''P{d+d-a + d+ad^a)} 
-I2{d+d-a + d+(yd-a){d±af 

+ { {-dip - 2id^pr) - \di + Idip + Id^pd^] 



32 



d^d^p 



.2p 



a + —b'{—ad^d^a + d^ad^a) 



(22) 



Usually the equation given by g^^ or g can be regarded as the constraint 
equation with respect to the initial or boundary values. The equations ob- 
tained here are, however, combinations of the constraint and a-equation of 
the motion since the tensor ^f^j, under consideration is the product of the 
original metric tensor in (|^) and the cr-function e"^"". 

Polyakov's anomaly induced action (|ll]) has the following form in the 
conformal gauge (]TP|) as usual 



r{o) 



127r 



J (fxpd+d^p 



and the variations with respect to g'^^ are given by 



<5^±± 487r 



e^" (9±p)2) +t±(x^)e- 



2p 



(23) 



(24) 



Here t^(x^) is a function determined by the boundary condition. The "po- 
tential" term T^'^^ in (|1^) has the following form in the conformal gauge 
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r(V) = JL[ rf2^g2p f8^-2pQ Q + 2) ( 1 - In 

967r J ^ V 



6/i2 



(25) 



and the variations with respect to g^^ are given by 



5^ 



±± 



xe 



(26) 



The variations with respect to p are given by 



s /r 



ind 



6p \ An 



5rw 
dp 



r 64 64 / 

6 J _!!^e-2''(9+9„p)V + — (ae-2pa+a 
L 3 3 



r 1 28 

+6' |-32(a+9_a)2e-2^ - -^e-2p^^^_^(^^^^_^) 
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-169+9- a 
32 



+y [-2d+d-ae-^pd+d-p + a+<9_ (a+a^ae-^^)} 

6" + hh + h')\ [iQd+d^ [e'^P{d+d-a + d+aO^a)] 



-48e-2p(a+<9_a + a+da^a)^ 
67r 



487r 



a+a_p 

e2^(8e-2''a+a_p + 2) |^1 - In 

8e-2^a+a_p + 2\ 



(27) 
(28) 



8e-2^a+a_p + 2' 

6/i2 



d+d^ In ■ 



127r 

+—e2''a+a_p In 



6/i2 y 

Se-'^Pd+d.p + 2 
6/x^ 



(29) 
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The variations with respect to a may be found as 



5 /r 



ind 



Sa \ 4:71 



6a 



+b' [32a+a_(e-2P5+a_(T) 
64 / 

+— {f^Pd+d-pd+d-a + d+d-iae-^fd+d-p 

y + 

+y {a+(e-2p5+9_p)a_a + d.{e~^Pd^d^p)d+a] 

+ [iQd+d-p + y9+a_(e-2''a+9_p) 



6" + -(6 + [2a+a_(8e-2''9+a_p + 2) 

-2 {a+(8e-^^9+(9_p + 2)d-a + a_(8e^2''a+a_p + 2)d+a] 
+489+(9_ {e-2''(a+a_a + 
-485+ {e-2p9_(T(a+a_a + a+aS^a 
-48a_ [e-'^Pd+a{d+d^a + a+aa_a 



5a 







(30) 
(31) 



It often happens that we can drop the terms hnear to a in ([T7|). In particular, 
one can redefine the corresponding source term as it is in the case of IR sector 
of 4D QG [ 



5 



±± 



L ind 



In that case, Eqs. (p2D, ( 127|) and (^) are reduced into 

+ -e^Pd±ad±{{8e~^Pd+d^p + 2)(r} 
3 

+^e^fad±{8e-^Pd+d^p + 2)d±a 



b" + -(b + 6') j [Ae^'d^ad^iSe-^Pd+d-p) 
A{d±afd+d-p + l2e^Pd±ad±{e~^P{d+d-a + a+aa_a)} 
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6p V 47r 



s /r, 



ind 



Sa V 47r 



-12((9+a_(j + a+a5_a)(a±cT)2' 

+ { {-dip - 2{d^pr) - \di + ^dip + ^d^pd^] 



X 



16 

—b'{—ad+d^a + d+ad^a) 



b" + -{b + b')^{d+d-a + d+ad.a^ , (32) 
b' I -32(a+a_ (7)^6-2'' - ^e-^Pd+d_p{ad+d_a) 

-— {-2d+ad-ae-^Pd+d-p + d+d- {d+ad.ae-^')} 
- |6" + ^(6 + 6')} [iGd+d- {e-^^id+d-a + d+ad^a)} 



-iSe-^^id+d-a + d+ad-af] 
b' [32d+d-{e-^''d+d-a) 
+y (e-^Pd+d.pd+d-a + d+d-iae'^^d+d-p] 

+^ {a+(e-2^a+a_p)a_(7 + d-{e-^fd+d-p)d+a} 



(33) 



2a+a_(8e-^''a+a_p+2) 

-2 {a+(8e-2^a+a_p + 2)a_a + d^{8e-^fd+d_p + 2)a+£7} 
+489+a_ {e-2''(a+a_(7 + d+ad_a)} 
-485+ {e-2''a_f7(5+a_(7 + d+ad-a)} 

-48a_ {e-2''a+(7(a+a_c7 + 9+(7a_(7)}] (34) 

The classical Einstein action 

1 



5' = 



IQttG 



rfSy^{i?(^)-2A} 



(35) 
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is reduced into 2 dimensional one in the conformal gauge (|T9|) as follows 
1 



S 
47r 



and we obtain 
6 / S 



(fx 



7±± 
6 



Sa \47r 



An 
S 



5_ /S_ 

6p \Att 



4e2"a+a_(p + a) + 2e^P+^'' - e^^+'^A - 4e2"a+a9_aJ , 

(36) 



16nG 
1 



(37) 



8e'" {39+a_a + 3a+aa_(T + d+d^p} 



levrG 
1 



4a+a_e'" - 2e'''+^"A + 4e='^+'" 



(38) 
(39) 



Above equations give the complete system of equations of motion for the 
system under discussion. 



4 Anti-evaporation of Nariai black holes 

Following the work by Bousso-Hawking we consider the Schwarzschild- 
de Sitter family of black holes, or more exactly its nearly degenerated case, 
so-called Nariai solution |Q. We are now interested in the instability of the 
Nariai solution. The Schwarzschild type black hole solution in the de Sitter 
space has two horizons, one is the usual event horizon and another is the 
cosmological horizon, which is proper one in the de Sitter space. The Nariai 
solution is given by a limit of the Schwarzschild-de Sitter black hole where 
two horizons coincides with each other. In the limit, the two horizons have 
the same temperature and they are in the thermal equilibrium. Near the 
limit, however, the temperature of the event horizon is higher than that 
of the cosmological one and we can expect that there would be a thermal 
flow from the event horizon to the cosmological one. Therefore the system 
would become instable and the black hole would evaporate. We also have to 
note that above cosmological black holes may naturally appear only through 
quantum pair creation which may occur in the inflationary universe . 
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In the Nariai limit, the space-time has the topology of x S"^ and the 
metric is given by 

ds^ = J (sin^ xdiJ^ - dx^ - dn) . (40) 

Here the coordinate x has a period vr. If we change the coordinates variables 
by 

r = In tan — 
2 



t 
4 ' 



(41) 



we obtain 



ds"^ = ^ ( -rft^ + dr^) + -dn . (42) 



A cosh^ X ^ ^ A 

This form would correspond to the conformal gauge in two dimensions. Note 
that the transformation (|4l|) has one to one correspondence between [i/j, x) 
and {t, r) if we restrict x by < x < ^i" (^^ runs from — oo to +oo). 

Now we solve the equations of motion. Since the Nariai solution is char- 
acterized by the constant 0(or cr), we now assume that a is a constant even 
when including quantum correction; 

cr = (To (constant) . (43) 

We also first consider static solutions and replace d± by ±1^^. Then the 
total constraint equation obtained by (^4|) , (p6|) , (|32D and (|37| ) becomes 



= -^e'^{dlp-{d,pf)+t,e'^ 

2„- -2e-2^a2p + 2 .... 
xe^^'ln . 44) 



The equations of motions given by p8|), (129|), (|33|) and (|39D (using the vari- 
ation of p) (0), (|34|) and (|38|) (using the variation of cr) have the following 



13 



forms: 



47rG V 



2 



247r V 6/i 



47rG' 

2 



+167r \^h" + -{b + 6')} 52 (e"2^a2p) . (46) 
Assuming a solution is given by a constant 2d scalar curvature 

R = -2e~^Pd'^p = Ro (constant) (47) 

and substituting ( ^7D into (^^) and (^B|), we find that Eqs. ( |45D and (^) can 
be satisfied if ctq and i?o satisfy the following two algebraic equations 

= Ro - AAe^"" + 4: (48) 







^(-Ae^- + 2e2-) + ^(-i?^ + 4)ao 



2nG V J 3 

247r " V y y ^ 

The above equation can be solved with respect to ctq and Rq in general 
although it is difficult to get the explicit expression. A special case is given 
by dropping the term corresponding to the terms linear to a in (|1^ (i.e., 
by dropping the second term proportional to b in (^9D ) and by choosing the 
renormalization scale /j, to be 

H = . (50) 



In this case, the solution of (|48|) and (|^) is given by 

I NGA 
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If we choose + sign in front of the root in the above equations, the solution 
reduces to the classical solution in the classical limit of — > 0. 
Equation (^Tj) can be integrated to be 

idrpY = -^e'^ + C . (52) 

Here C is a constant of integration. Substituting p7|) and (|5^) into (|4^), we 
find Eq. ( PD is satisfied if and only if 

C. (53) 



19277 

Eq.(|52|) can be integrated to become 



Ro cosh' (rVC^' 



We now consider the perturbation around the Nariai type solution (^) 
and (H); 

p = Pq + eR{t, r) 

to = -^C + eT^ . (55) 
1927r ^ ' 

Here e is a infinitesimaly small parameter and po is a solution given by (|5^): 

g2po ^ ^ ^ C55) 

Rq coshVrv^V 



Then we obtain the following equations 
cosh^ {r^/cY^^ 



= -Stt -^—rdiS 
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^4.(-f.V„-|(. + 6')) (j 



Cd.d^S 



=^ (~R + dlR ± tanh (ry/c) d±R 
2o V 2 =^ ^ V 



TV C 



247ri?n 



2C 


1 


fNC 


-^0 cosh^ 









N 



C 



dlR- 



In 



i?o + 2 



+ 



247r |_2i?o ^ cosh^ (r^/c) 6//2 



-Ro cosh^ (rv^ 
8C f /S 1 



d^d-R 



Ro\\^ cosh" (r^fc) 2 

T tanh (rVc) d±d+d_R^ 



1 2C 



IGttG i?o cosh^ (rVCj ^ \ J J 

= 47r6'|-yi?o(To9+a_5+y^a+a_(cosh2(r^)5+9_S)| 

-47r(6 + (cosh^ (r^) 

AT C 1 



127r i?o cosh^ (r VC) 
Ro + 2' 



(57) 



X <^ 2i? 1 - In 



+ 



127ri?o + 2 I 
1 



-2Rod+d_R + 
8d.d_S - 



Ro + 2 



C 

SAC 1 



Ro cosh^ (rVc) 



2RoR + ^ cosh^ (rVc) 5+5- | 
d+d_ (cosh^ (rVc) 

e2"°(i? + 25) 



16 



16C 



(R + S) 



Ajrb' S^^Rod+d^S + 16^9+a_ (cosh^ (r^C) d+d^S 

d+d- (cosh^ (r^/c) d+d^R 
-47r X 16(6 + b')^d+d- (cosh^ (rVc) d+d^S 



(5^ 



16 



jaoRod+d^R + Y~c~ 



1 



WttG 
16AC 



8 3d+d_S + d+d_R + 



C 



Ro cosh^(rv^) 



{R + 2S) + 



2 cosh^ (r^c) 
16C7 1 



S 



Ro cosh^ (ry/c) 



(R + S) 



(59) 



When deriving the above equations, we have dropped the terms hnear to a 
in (0), for simphcity. 

The equations (|58|) and (|59|) can be solved by assuming i? and S are given 

by 



/2(t,r) = Pe^*°^cosh"ryC 
S{t,r) = Qe^'"^ cosh'' rVC . 



(60) 



Here P, Q and a are constants. We should note that a < if we require 
R and S are finite in the limit of r ±oo. Then we obtain the following 
equations 



cosh^ (^rVc) = AR 

cosh^ (ry/C)d+d^S = AS 

a{a — 1)C 



A 



(61) 



Note that there is one to one correspondence between A and a if we restrict 
^ > and a < 0. Using (|6lD, we can rewrite (^) and ( ^9]) (using (|48|)) as 
follows: 







3 [~ ^ ~c] 3 ~C~~AtiG' 



e2-o ^2A -C))Q 
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f N C I, , -R„ + 2\ N C / 2 ; 



QA-C - — ])Q 



In order that the above two algebraic equations have a non trivial solution 
for P and Q, A should satisfy the following equation 

f N C , i?o + 2\ iV C7 / 2 ; 

+^^r-A+4A^i-^e-r-i+A 



67r /2o + 2 V C / SttG V i?o 
X S^GiTTb'Ro (^A + 1^2) - 64n{b + b')^A^ 



'''' ^,A-C-'^ 



AtxG V Rq . ^ 

^ F{A). ' (63) 



If the equation ( pSf ) has a positive (A > 0) solution for A, i.e., a real negative 
solution for a, there is a solution where S increases exponentially in time and 
the system would be unstable. It is generally difficult to find the parameter 
region where A has a positive solution but it might not be so difficult to show 
the existence of such a parameter region. 

In order to show the existence, we now consider the solution of (^Tj) and 
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Ro^i limit (A^GA 0). Then F(0) is given by 



m - ^ > . (64) 



c 



On the other hand, when A — ^ 



^ 2 / 247r V 23 ■ 32 ■ Stt 27rAG / ^ ' 



Eqs.(|64D and (|65[) tell that there is a solution of the equation -F(A) = when 
< A < f , i.e., 

> a > -1 . (66) 

In order to show that the existence of the solution (^) implies the instability 
of the system, we consider a linear combination of (|60|): 



S{t, r) = Q cosh taVC cosh" rVC . (67) 

It should be noted that any linear combination of two solutions is always 
a solution since the perturbative equations of motion are linear differential 
equations. The horizon is given by the condition 

Va • Va = . (68) 

Substituting (|67|) into (|68|) , we find the horizon is given by 

r = at. (69) 

Therefore on the horizon, we obtain 

S{t, r(t)) = Q cosh^+" taVC . (70) 

It should be noted 1 + a > in the solution of (§§), which implies the system 
is unstable. Since the radius of the horison Vh is given by 

rh = e'' = e-^o+^^^^'^W) , (71) 

the radius increases monotonically in time if the initial perturbation at t = 
is positive Q > 0, on the other hand, if the initial perturbation is negative 
Q < 0, the radius shrinks, i.e., the black hole evaporates. 
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The solution in > a > — 1 shows that the system is unstable. The 
perturbation, however, becomes stable if there is a solution where a < — 1, 
i.e., A > ^. In order to show the existence of the solution where A > ^, we 
consider F{A) in the limit of ^4 ^ +oo. As in the previous case, we consider 
the limit near the classical solution Rq ~ 4. Then we find 

Combining (^) and ([7^), we find that there is a solution where A > ^ i.e., 
a < — 1. This implies that the perturbation is stable since the perturbation 
becomes exponentially small in time on the horizon ([fO|) 

5(t,r(t)) ^Qe(^+°)*l"l^ . (73) 

Therefore if the initial perturbation is negative Q < and the perturbed 
radius of horizon is smaller than that of the Nariai limit, which is the thermal 
equilibrium, the radius increases in time and approaches to the Nariai limit 
asymptotically, i.e., the black hole anti-evaporates as observed by Bousso 
and Hawking p[. It might be interesting that there are both of stable and 
unstable perturbations in the model discussed here although the treatment of 
the quantum corections is rather different from that of Bousso and Hawking. 
In the model by Bousso and Hawking, it depends on the initial conditions 
if the perturbation is stable or unstable, i.e., the perturbation is stable and 
the black hole anti-evaporates if its initial value does not vanish but its time 
derivative vanishes, on the other hand, the perturbation is unstable if its 
initial value vanishes but its time derivative does not vanish. In the model 
discussed in this paper, however, the initial conditions for the value and the 
time derivative of the stable and unstable perturbations are completely the 
same, i.e., the initial value does not vanish but the initial time derivative 
vanishes for both cases. The reason why these initial conditions lead to the 
physically different results could be because the effective Lagrangian contains 
higher derivatives and the initial value and the initial time derivative are not 
sufficient as a set of the initial conditions and we need the initial conditions 
for the higher derivatives to determine the dynamics. 
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5 Discussion 



In summary, we demonstrated that in the model under discussion black holes 
may evaporate or anti-evaporate. We have to note that our consideration has 
quite common character. We limited to the case of conformal scalars but we 
could start from the theory of free conformal scalars, spinors and vectors. 
(As is known many asymptotically free GUTs in the early Universe are 
also asymptotically conformally invariant and they may be represented as 
collection of free conformal fields). Then, the anomaly induced effective 
action W would have the same form as in section 2 where only coefficients b 
and b' would be changed (but their signs are the same). And in the Eqs.([TI|) 
and ([T^ ) again only numerical coefficients would be changed due to spinor 
and vector corrections (in Eq. ([T^) also signs of some coefficients could be 
changed). Hence, our result actually corresponds to quantum conformal 
scalar-spinor-vector system. 

There are different possibilities to extend our results. From one side, one 
can study the similar, anti-evaporation effect in black holes realised in gravity 
with torsion. The corresponding anomaly induced effective action is given 
in last ref. of 0. That could help to find finally some phenomenon where 
torsion effects may be significally strong to prove (or disprove) the existence 
of torsion in the Nature. 

From another side, using the effective action of section 2 one can inves- 
tigate the quantum cosmology of Kantowski-Sacks form [|1^. Recently, such 



cosmology has been studied in large and s- waves approximation [0. In 
the above works it was found that realization of non-singular Universe for 
the models under discussion is problematic. The results of the present work 
(if exchange time and radius) indicate to the possibility of construction of 
non-singular Kantowski-Sacks early Universe. 
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